Two-state theory of nonlinear Stochastic Resonance 
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An amenable, analytical two-state description of the nonlinear population dynamics of a noisy 
bistable system driven by a rectangular subthreshold signal is put forward. Explicit expressions 
for the driven population dynamics, the correlation function (its coherent and incoherent part), the 
signal-to-noise ratio (SNR) and the Stochastic Resonance (SR) gain are obtained. Within a suitably 
chosen range of parameter values this reduced description yields anomalous SR-gains exceeding unity 
and, simultaneously, gives rise to a non-monotonic behavior of the SNR vs. the noise strength. The 
analytical results agree well with those obtained from numerical solutions of the Langevin equation. 



The phenomenon of Stochastic Resonance (SR) at- 
tracts ever growing interest due to its multi-facetted rel- 
evance for a variety of noise-induced features in physics, 
chemistry, and the life sciences 0, 0, El S 01 • Several SR- 
quantifiers have been used to characterize the response 
of a noisy system to the action of time-periodic external 
forces. In particular, the non-monotonic behavior of the 
output signal-to-noise ratio (SNR) with the strength of 
the noise has been used widely. A dimcnsionless quantity 
that measures the "quality" of the response with respect 
to the input signal is the SR-gain defined as the ratio of 
the output SNR over the input SNR. Ideally, one would 
wish to obtain the characteristic amplification of the SR 
phenomenon 0, El and, simultaneously, SR-gains larger 
than unity. For superthreshold sinusoidal input signals 
SR-gains larger than unity have been reported before Q • 
In recent analog 0, EJ [l^ and numerical 0, 0] simu- 
lations of noisy bistable systems driven by subthreshold 
multifrequency input forces, surprisingly large SR-gains 
larger than unity have been established. 

In order to clarify the conditions under which these 
anomalous large SR-gains occur, it would be interest- 
ing to propose simplified models, amenable to analytical 
treatment, which describe this rich behavior of the re- 
sponse. A detailed proof that SR-gains larger than unity 
are incompatible with Linear Response Theory (LRT) 
has been presented in Thus, any theoretical expla- 
nation of the simultaneous existence of SR and anoma- 
lous large gains is rooted in the response beyond LRT. 

The main focus of this work is to present such a sim- 
plified two-state description of the nonlinear dynamics of 
a noisy, symmetric bistable system driven by a rectan- 
gular subthreshold signal. A two-state description of SR 
has been considered previously in the pioneering work 
by McNamara and Wiesenfield for sinusoidal input sig- 
nals |l3|. In clear contrast to their work, however, we 
will here not linearize the transition probabilities in the 
strength of the applied force. In doing so, we put forward 
explicit analytical expressions for the driven population 
probabilities, the nonlinear correlation function (its co- 



herent and incoherent part), the SNR and the SR-gain. 
These novel nonlinear response results come forth solely 
because the rectangular signal - in contrast to sinusoidal 
driving - involves two force values only. 

To start, let us consider a system characterized by a 
single degree of freedom, x, whose dynamics (in dimen- 
sionless units) is governed by the Langevin equation 



x(t)=-U' [x(t),t]+t{t), 



(1) 



where £(t) is Gaussian white noise of zero mean with 
(£(t)£( s )) = 2DS(t—s), and — U'(x, t) represents the force 
stemming from the time-dependent, archetype bistable 
potential 



U(x,t) = ^-^-F(t)x. 



(2) 



We will restrict our study to a periodic rectangular driv- 
ing force with period T, 



F(t) = (-l) n(t) A, 



(3) 



where n(t) = [2t/T\, [z\ being the floor function of z, 
i.e., the greatest integer less than or equal to z. In other 
words, F(t) = A [F(t) = -A] if t e [nT/2, (n + l)T/2) 
with n even (odd). Our focus is on subthreshold signals; 
more precisely, we will assume that A < A t h — v/4/27 
where A t h is the static threshold value (the dynami- 
cal threshold value always exceeds this adiabatic thresh- 
old A t h)- In this case, the potential U(x,t) presents 
two minima at x%(t) < and X2(i) > 0, and a maxi- 
mum at XM{t)- Because the potential fulfills the sym- 
metry property U(x,t + T/2) = U(—x,t), then XM{t) = 
(-l) n <*W (0), and 



Xj(t) 



-1) 



, MQ) 



(-i) 



xm(0) 



(4) 



where Ax(0) = £2(0) — .xi(0). Additionally, we have 
taken into account the symmetry relation x\ (t) + x^it) + 
x M {t)=0. 
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Throughout the following, we will assume that the 
noise strength, D, is sufficiently small so that the in- 
trawell relaxation time scale is negligible compared with 
the time scale associated to the interwell transitions and, 
as well, the driving time scale T/2. In this case, and for 
sufficiently slow driving 0, 0] , the exact Langevin dy- 
namics can be approximated by a nonstationary, Marko- 
vian two-state description of the form 

P(l,f) = - ll (t)P(l,t)+ l2 (t)P(2,t), (5) 
P(2,f) = -7 2 (t)P(2,i)+ 7l (*)P(M). (6) 

Here, P(l,i) and P(2,t) denote, respectively, the pop- 
ulations to the left and to the right of the maximum 
position, XM(t), and 7j(f) is the Kramers' rate of escape 
|l6| from the well j at time t, i.e., 



7i(*) 2n exp 



U[x M (t),t]-U[ Xj {t),t] 
D 



j , (7) 

where LJj(t) = y/U"[ Xj {t),t] = y/3[xj(t)] 2 - 1 and 

w M (i) = VP"[xM(t),t}\ = ^/l-3[x M (t)} 2 = ^m(O). 
According to the above mentioned symmetry property of 
the potential, 7 j(f) can be expressed in the form 



7i(*) = ^ 



1 _ (-l)"( t )+^ AP e9 (0) 



(8) 



where T = 71 (0) +72(0) = 71(f) + 72(f), and AP eq (0) = 
P eq (2, 0) - P eq {l, 0), Peq{j, 0) being the equilibrium pop- 
ulation of the state j corresponding to the rates taken 
at time t = 0, i.e., P eq (J, 0) = [5 jA 72 (0) + 5 jt 2 71 (0)] /T. 
Eqs. © and © describe the evolution of the popula- 
tion between two consecutive changes of shape of the 
potential. These populations can be discontinuous at 
t = nT/2, as a result of the sudden change in the location 
of the maximum at those instants of time. Therefore, in 
order to complete the description, jump conditions for 
the populations at t = nT/2 have to be added. For our 
situation with a rectangular signal we shall assume that 
the probability distribution before a change of the po- 
tential is sufficiently well localized around the minima 
that there is almost no probability transfer from one well 
to the other as the maximum location changes. In this 
case, the populations can be considered to be continuous 
at t = nT/2, i.e., 

km P(j,f) = lim P(j,t)=P(j,nT/2). (9) 

This relation becomes exact at asymptotic weak noise. 
At finite weak noise strength D those finite jump con- 
ditions induce for the correlation function some small, 
weakly nonanalytic structures which in turn may cause 
dip-like features around even numbered multiples of the 
angular driving frequency £1. Indeed, such dips in the 
spectrum are well known to occur at weak noise for time- 
continuous, e.g. sinusoidally rocked bistable systems in 



the nonlinear response regime |l7|: note, however, that 
for our case with a constant force (or potential) the source 
of the nonanalytic structure is of different origin. For 
rectangular driving forces and noise strength values lead- 
ing to SR-gains larger than unity, the incoherent part 
of the correlation decays very rapidly on the time scale 
T/2 so that these small corrections stemming from the 
approximation in Q can safely be neglected. After us- 
ing the normalization condition Pi(l, f) + Pi (2, f) = 1, 
Eqs. and ijfj} yield 



p(i,f) = -rp(i,f)+ 72 (f). 



(10) 



Upon observing that 72(f) remains constant between 
two consecutive changes of the potential, the solution of 
Eq. (fTUf) can be expressed as 



P(l,f) 



72(f) 



P[\,n(t)T/2] 



72(f) 

r 



-r t 



(11) 

With this relation, it is straightforward to evaluate the 
time-periodic, asymptotic long-time solution of Eq. llUfl , 
i.e. P 00 (l,f). In order to do so, we calculate the values 
P oc (l,nT/2) by making use of the symmetry property 
P oc (l,f ±T/2) = Poo(2,f) = 1 -P 00 (l,f), as well as of 
Eq. After inserting the result in Eq. (|11|) . one finds 
the T-periodic solution 

^(1, f) = ~ [l - (-l)" (t) AP e(? (0) 

-r[t-sipL] 

+(-l)«W AP eg (0) — . (12) 

1 + e 2 

The average of the coordinate in the long-time limit, 



(x(t)) 



(TS) 



Si=i x j{t)Poo(ji f)i is evaluated to read 



Mf)>£ 5) = (-ir^^(x(0)) eq -[x M (0) + 2(x(0)) ( 
-r[t-^] ^ 



1 + e~T- 



(13) 



where (x(0)) eq = £'=i Xj (0)P eq (j,0). 

The conditional probability P(l, f \j, to) can also be cal- 
culated using the observation that P(l,t\j, f ) — Poo(l,f) 
(for f > f ) fulfills the homogeneous equation obtained 
by removing the term 72(f) from Eq. (|10f) . with the ini- 
tial condition P(l, t Q \j, f ) - Poo(l, f ) = <S,-,i -P oo (l,f ). 
Thus, the result for f > f is 

P(l,f|j,f ) =Poo(l,f) + [^,i-Poo(l, fo^e^'-* ), 

(14) 

and, likewise, P(2, t\j, t ) = 1 — P(l,f|j, f ). In terms 
of the time-periodic one-time probability in Eq. (|12l) the 
two-time joint probability reads 

Poo(j,f;Mo) = (-l) J+fe Poo(l,fo)/ 3 oo(2,fo)e- r (*- t °) 
+P oo (j,f)P oo (fc,f ), (15) 
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for t > to. Therefore, the asymptotic two- 
time correlation function, (x(t)x(to))aa = 



LLi Efe=i x i CO 2 * (*o) -Poo (j, t;k,t ), reads 



(x(i)a;(to))^ 5) = [Ax(0)] 2 P 00 (l,io)Poo(2,io)e- r ( t -*«) 
+<z(*)>^<*(i )}£ S) , (16) 

for t > to, where we have used that Ax(t) = Ax(0). 

The two-time correlation function (x(to + i")x(io))oo 
is a periodic function of to with the period of the exter- 
nal driving 6]. Then, it is convenient to apply a time- 
average to obtain the time-homogenous correlation func- 
tion, C< ts )(t), i.e. 

C(TS)(T) = f£ dt o(^o+r)x(t ))^ s \ (17) 

In virtue of Eq. (JS), C (ts \t) can be written as the 
sum of two contributions: a coherent part, C^^\t), 
which is periodic in r with period T, and an incoher- 
ent part, C^J h (r), which decays to as r — > oo. These 
two contributions are given by, respectively, 

Cg?(r) = ± jfd*, (x(t + T))g s Hx(t ))g s \ (18) 



where fl = 2ir/T is the angular frequency of the external 
driving, and H(u>) denotes the Fourier cosine transform 
of H(t), i.e., H(u) = 2/ir J™ dr H (r) cos(wr). Note that 
this definition of the output SNR differs by a factor 2, 
stemming from the same contribution at lo = —Q, from 
the definitions used in earlier works 0, @ • The period- 
icity of the coherent part gives rise to delta peaks in the 
spectrum. Thus, the only contribution to the numerator 
in Eq. I)22ll stems from the coherent part of the correla- 
tion function. The output SNR can then be expressed 



R 



(TS) _ Qu 

Hi 



(23) 



where 



and 



Q { J S) = J / drCfJ\r)co S (nr), ,24) 



f S) = zl drC^^cos^r). (25, 



Then, from Eqs. $ZQ, JUJ, and we get after- 

some simplifications 



and 



( TS) , , = [A^(0)] 2 e 



-IV 



dt Poo(l,*o)-Poo(2,t ). 



(19) 

Upon combining Eqs. lJT2)l. (JT2J), ((El and (HHJ, one 

obtains after some cumbersome algebra 



C 



(TS) 



4r 

2»(t) + 1 - — 



h M = (-i) nw |(^(o))^ 

4sech (rf ) sinh { ^ [2n(r) + 1 - *f ] } 



and 



incoh\ / 



xm(0) 



[Ax(0)} 



TT 



(20) 



2 „-r T 



4tanh(^) 



rr 



l+[AP eg (0)r 



- 1 



(21) 



According to McNamara and Wiesenfeld 01 > the output 

(TS) 

SNR, R out , is defined in terms of the Fourier transform 
of the coherent and incoherent parts of (r) as 



n(TS) 



(TS) 



(22) 



Q 



(TS) _ 



2 2 



4(x(0)) 2 g r 2 + [x M (0)]^ 



and 



,(TS) 



7r 2 (r 2 + f! 2 ) 



fAx(o)l 2 r rA , Sl2 

1 > 1 + [AP eq (0)j 2 



(26) 



2tt (r 2 + n 2 ) 



4tanh(^) 

rr 



(27) 



The signal-to-noise ratio of the input signal, F(t) + £(t), 
can readily be evaluated from the definition, yielding 



R ~ 4A 
ttJJ 



(28) 



Thus, the SR-gain which is defined as the ratio of the 
SNR of the output over the SNR of the input, emerges 
as 



G (TS) = 



R. 



(TS) 



R 



(29) 



nip 



and it can be evaluated explicitly upon combining the set 
of Eqs. (E3J), (H3), and JUJ. 

In Fig. 1 we compare our analytical results for the be- 
havior of several SR-quantifiers as a function of the noise 
strength D with results obtained by numerically integrat- 
ing the Langevin equation [Eq. l|T|l]. The numerical pre- 
cise solution is based on an algorithm due to Greenside 
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FIG. 1: Stochastic Resonance beyond linear response: Sev- 
eral SR-quantifiers are depicted vs. the noise strength D; 
namely, the cosine transform of the coherent part of the driven 
correlation function at angular driving frequency SI, i.e. nu- 
merator of the SNR (Q u ), and the cosine transform of the 
incoherent part, respectively, i.e. the denominator (Qi), the 
output SNR (Rout) and the SR-gain (G). These character- 
istic quantities are evaluated for a rectangular driving force 
with angular frequency SI = 0.01 and subthreshold amplitude 
A = 0.25. The solid lines depict the results obtained within 
the two-state description, whereas the numerical precise re- 
sults obtained from the driven Langevin dynamics in 11I2H are 
given by the circles. 



and Helfand; for details see in the Appendix of Ref. 01 • 
The value of the subthreshold input amplitude is chosen 
as A = 0.25, and the frequency of the external driving is 
set at ft — 0.01. This small value of the angular frequency 
has been chosen in order to observe a characteristic non- 
monotonic behavior versus D of several quantifiers associ- 
ated to SR and, simultaneously, SR-gains exceeding unity 
over a wide range of D (see Ref. |l2|). The agreement 
between the analytical nonlinear response results and the 
numerical results is surprisingly good at moderate noise 
values, and becomes even excellent for small values of 
D. The better agreement at low noise values corrobo- 
rates with the fact that with increasing noise st reng th 
the Markovian two-state description also worsens [15j . 

With this work we have put forward an analytical two- 
state description for both, the nontrivial, nonlinear pop- 
ulation dynamics and the nonstationary correlation be- 
havior of noisy bistable systems driven by periodic rect- 
angular subthreshold forces. It is indeed remarkable that 
our analytical, nonlinear two-state approach does capture 
well, both the non-monotonic, bell-shaped behavior of 
the nonlinear SNR vs. noise strength D, i.e. the charac- 
teristic SR phenomenon, and the occurrence of SR-gains 
larger than unity. This latter result is a true benchmark 



of the nonlinear response behavior of a driven bistable 
stochastic dynamics. Our analytical findings corroborate 
those obtained by means of numerical solutions of the 
Langevin equation [l2^ . 

This very two-state theory beyond linear response of a 
driven, metastable stochastic population dynamics likely 
proves useful also for phenomena other than Stochastic 
Resonance: it equally well can be applied and generalized 
to describe the behavior of rocked Brownian motors |l8| , 
the description of the switching dynamics over adiabati- 
cally sloshing potential landscapes or also to driven 
noisy threshold characteristics in general, such as e.g. for 
a driven neuronal noisy spiking dynamics. 

The authors acknowledge the support of the Direction 
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03822), the Junta de Andalutia, DAAD-Acciones In- 
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